Electric Dipolar Kondo Effect Emerging from Vibrating Magnetic Ion 
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When a magnetic ion vibrates in a metal, it inevitably introduces a new channel of hybridization 
with conduction electrons and in general, the vibrating ion induces electric dipole moment. In 
such a situation, we find that magnetic and non-magnetic Kondo effects alternatively occur due to 
the screening of spin moment and electric dipole moment of vibrating ion. In particular, electric 
dipolar two-channel Kondo effect is found to occur for weak Coulomb interaction. We also show that 
magnetically robust heavy-electron state appears near the fixed point of electric dipolar two-channel 
Kondo effect. We believe that the vibrating magnetic ion opens a new door in the Kondo physics. 
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It has been widely recognized that Kondo phenomena 
generally appear when localized entity with internal de- 
grees of freedom is coupled with conduction electrons. 
Concerning the original problem of resistance minimum 
phenomenon in metals with magnetic impurities, Kondo 
has actually shown it by quantum-mechanical calcula- 
tions for scattering amplitude of electrons due to mag- 
netic impurities [l|. Then, it has been revealed that the 
singlet state is formed from local magnetic moment due 
to the coupling with conduction electrons y. After the 
understanding of the Kondo effect in the dilute magnetic 
impurity system, interests of researchers have moved to 
the impurity with complex degrees of freedom. 

One research direction has been found in the explicit 
consideration of orbital degree of freedom of localized 
electron. Coqblin and SchriefFer have derived exchange 
interactions from the multiorbital Anderson model [3|. 
Then, the concept of multi-channel Kondo effect has been 
developed on the basis of such exchange interactions [J] , 
as a potential source of non-Fermi liquid phenomena. 
Such non-Fermi liquid properties have been pointed out 
also in a two-impurity Kondo system [5|, |6[ . Concerning 
the reality of two-channel Kondo effect, Cox has pointed 
out the existence of two screening channels in the case of 
quadrupole degree of freedom in a cubic uranium com- 
pound with non-Kramers doublet ground state |7| . 

Another possibility is non-magnetic Kondo effect with 
phonon origin. First Kondo has considered Kondo-like 
behavior in a two-level system [8|. The two-level Kondo 
system has been proven to exhibit the same behavior as 
the magnetic Kondo effect [9[ . Yu and Anderson have dis- 
cussed the phononic Kondo effect from a different view- 
point [10|. They have pointed out that the scattering 
process between spinless s-wave conduction electron to 
p-wave one is produced by ion displacement. 

Recently, the Kondo effect with phonon origin has at- 
tracted renewed attention due to active experimental in- 
vestigations on cage structure materials, in which a guest 
ion is vibrating in a cage composed of relatively light 
atoms. We believe that the vibration of magnetic ion in 



the cage provides a new ingredient in the Kondo physics, 
since dynamical aspects of magnetic impurity have been 
considered unsatisfactorily in the Kondo problem. In 
fact, quite recently, the two-channel Kondo effect has 
been confirmed in the model for vibrating magnetic im- 
purity lll-ll4| . Also for the promotion of our understand- 
ing on magnetically robust heavy-electron phenomenon 
observed in cage compound [l5|, we further develop the 
Kondo physics of vibrating magnetic impurity. 

In this Letter, we analyze a two-channel conduction 
electron system hybridized with vibrating magnetic ion 
by using a numerical renormalization group technique. 
We confirm magnetic and non-magnetic Kondo effects 
originating from the screening of spin and electric dipo- 
lar moments, respectively, by evaluating entropy and sus- 
ceptibilities for spin and electric dipole moments. Near 
the fixed point for electric dipolar two-channel Kondo 
effect, we find magnetically robust heavy-electron state 
from the direct evaluation of the Sommerfeld constant. 

Let us consider a two-channel conduction electron sys - 
tem hybridized with vibrating magnetic impurity |12h14| . 
In the unit of 7i=fcB = l; the Hamiltonian is given by 



H 



E h(4..^ + h-c-) + 9x{cl,J, + h.c.)] (1) 



+ Un^ni^ + Efin^ + ni)+Lox'^/2+p'^/2, 

where Sk denotes conduction electron dispersion, Ckicr in- 
dicates the annihilation operator for conduction electron 
with momentum fc, angular momentum £, and spin a, 
fa^ is the annihilation operator for localized electron with 
spin CT, na=flfa, U is the Coulomb interaction between 
localized electrons, Ef denotes the local /-level energy, 
Vq is the hybridization between s-channel conduction and 
localized / electrons, g is the electron- vibration coupling, 
X denotes the ion displacement, p indicates the corre- 
sponding canonical momentum, and w is the vibration 
frequency. Note that we set the reduced mass of vibra- 




FIG. 1: (Color online) Schematic view of the situation de- 
scribed by the model Hamiltonian (1). 



tion as unity. In order to consider the symmetric case, 
we set Ef~—U/2 throughout this paper. 

In this paper, we analyze the model with the use of 
a numerical renormalization group (NRG) method J16l |. 
We introduce a cut-off A for the logarithmic discretiza- 
tion of the conduction band. Due to the limitation of 
computer resources, we keep N low-energy states. In 
this paper, we set A==5 and iV=:5000-- 10000. Note that 
the temperature T is defined as r=A^(*~^)/2 in the NRG 
calculation, where i is the number of the renormalization 
step. The phonon basis is truncated at a finite number 
A^ph, which is set as A^ph=300 in this paper. 

It is convenient to introduce phonon operators a and 
a^ through the relation of x^{a + a^)/\/2uJ. Then, we 
define the magnitude of phonon-assistcd hybridization as 
Vi=g/\/2uJ. The energy unit is the half of the conduction 
electron band Z), which is set as unity. As for parameters, 
we set Vb=w=0.2 and we change the values of U and Vi. 

First let us visualize the situation described by the 
present model. In Fig. 1, the two horizontal lines sym- 
bolically denote s and p channels which are hybridized 
with magnetic impurity. As shown in iJ, s-channel elec- 
trons are hybridized with localized electron in a stan- 
dard manner, while the hybridization process between 
p-channel and localized electrons is assisted by phonons. 
Note that parities of s- and p-channel electrons are dif- 
ferent, as easily understood from the values of angular 
momenta. Namely, s- and p-channel electrons possess 
even and odd parities, respectively. Since the parity for 
ion displacement x is odd, it is natural that the phonon- 
assisted hybridization occurs only for p channel. 

When magnetic ion stops at an origin, we consider only 
the screening of impurity spin moment by s-channel elec- 
trons, leading to the conventional Kondo effect. How- 
ever, when the ion is vibrating as shown in Fig. 1, there 
occurs another screening process due to p-channel elec- 
trons. In such a situation, in addition to the screening 
of spin moment, the electric dipole moment in propor- 
tion to the displacement x should be also screened by 
conduction electrons, leading to non-magnetic Kondo ef- 
fect. Intuitively, we understand that the main screening 
channel is converted between s and p due to the balance 
between Vq and Vi. In fact, the phases with different 
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FIG. 2: (Color online) (a) Phase diagram on the Vi-U plane 
for Vo=uj=G.2. The color gradation indicates the magnitude 
of average displacement q. Solid curve with solid circles de- 
note the two-channel Kondo line. Dotted curve with open tri- 
angles and broken curve with open squares denote d?q/dU^—Q 
and d^q/dVi=0, respectively. As for the definition of q, see 
the main text, (b) Entropy vs. temperature. 



quantum numbers in relation with parity have been found 
to be converted between the regions of large and small 
Vi. Then, two-channel Kondo effect has been confirmed 
to occur just at the boundary between those two phases 
11|-|13|. However, another conversion between magnetic 



and non-magnetic Kondo effects controlled by the bal- 
ance between U and Vi has not been clarified yet. Then, 
in this paper, we unveil such a new point in relation with 
magnetically robust heavy electron state. 

Now we explain the phase diagram in Fig. 2(a). As 
mentioned above, the solid curve indicates the bound- 
ary between two phases with different quantum numbers, 
which are obtained by the s- and p-channel Kondo screen- 
ing, respectively. On the boundary curve between two 
phases, the two-channel Kondo effect is realized. Since no 
phase conversion occurs on the line of C/=0, the boundary 
curve asymptotically approaches the line of U~0. Note 
that low-energy spectra of three fixed points of s-channel, 
p-channel, and two-channel Kondo states have been re- 
vealed in Refs. 12| and [13|. New results of this paper 
are color gradation and other broken and dotted curves. 
Their meanings will be discussed later. 

Next we discuss the NRG results of entropy. In 
Fig. 2(b), we show the change of entropy along the lines 
of C/=2, Vi=0.2, and t/=0.05. On the line of C/=2, for 
1^1=0.14 and 0.2, we find the plateaus of log 2, which 
are eventually released at low enough temperatures. At 
(Vi, C/)=(0.175,2), we observe the entropy of 0.51og2 
at low temperatures, which is the signal of two-channel 
Kondo effect. When we decrease the value of U from 
U=2 on the line of Vl=0.2, we observe the increase of 
the Kondo temperature Tk, which is characterized by 
the release of entropy log 2. This is quite natural from 
the viewpoint of the conventional magnetic Kondo effect. 
Thus, we deduce that the Kondo effect on the line of U=2 
as well as in the region of U^l on the line of Vi=0.2 is 
originating from the screening of impurity spin moment. 



However, for U^l, a plateau of log 2 again appears. The 
temperature region of log 2 is wider for smaller U and 
Tk is decreased with the decrease of U. Such behavior 
is contradictory to the conventional magnetic Kondo ef- 
fect. Furthermore, at (Fi, C/)=(0. 17384, 0.05), we again 
observe a clear plateau of entropy of 0.5 log 2, but it is dif- 
ficult to understand the origin of this two-channel Kondo 
behavior only from the result of entropy. 

In order to clarify what quantity is screened, we eval- 
uate susceptibilities for magnetic and electric dipole mo- 
ments, which are, respectively, defined by 
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dT{P{T)P). (2) 



denotes the operation to take thermal average, 
e"''Me-"^, M=gsHBini-ni)/2, and P=Ze{a + 



a')/\/2uj, where gs is the electron g-factor which is set as 
gs~2, ^B is the Bohr magneton, Z is valence number of 
ion, and e is electric charge. In the actual calculations, 
we normalize them as Xm/ij-b ^^'^ Xp/(-^^6^/2'^)- 

In Fig. 3(a), we show Txm for the same values of U 
and Vi in Fig. 2(b). As we have deduced above, on the 
line of U=2, we observe the decrease of Txm around Tk- 
Note that the curves of Txm for Fi=0.14 and Fi=0.2 at 
U^2 are quite similar at low temperatures, when T is 
rescaled by Tk- On the other hand, the curve of Txm 
for (Fl, C/)=(0.175, 2) is apparently different from those 
for Vi=0.14 and Vi=0.2. It is due to the non-Fermi hq- 
uid behavior in the two-channel Kondo effect, leading 
to xm^— logT. This logarithmic correction in xm is 
clearly observed. Along the lines of 1^=0.2 and f7=0.05, 
Txm decreases at relatively high temperature, but the 
release of the entropy does not seem to correspond to the 
temperature at which Txm is decreased. 

Now we turn our attention to the results for Txp in 
Fig. 3(b). The entropy release for (Fi, [/)==(0.2,0.05) 
corresponds to the decrease of Txp, indicating the oc- 
currence of the Kondo effect concerning electric dipole 
moment. Note that Xp is related to phonon Green's func- 
tion. In the strong electron-phonon coupling region, the 
center of oscillation is shifted either right or left, leading 
to xp(xq^/T, where q=\J{{a + at)2). Thus, for small U 
and large Vi region, we expect that Txp becomes con- 
stant at high temperatures, as found in Fig. 3(b). 

When T is decreased, Txp is decreased from the con- 
stant value due to the Kondo screening of electric dipole 
moment and it eventually goes to zero at T=0. This can 
be called the parity Kondo effect, since near degeneracy 
with different phonon parities characterizes the electric 
dipole, which is coupled with conduction electron parity, 
leading to the non-degenerate ground state with fixed 
total parity. Note that the total parity is specified by 
or 1, depending on U and Vi- At low enough temper- 
atures, we find Txp=2T /uji, where uji is renormalized 
phonon energy smaller than w. From the local Fermi- 
liquid theory, we find Tk cx; wi, but the proportional 
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FIG. 3: (Color online) Susceptibilities for (a) magnetic and 
(b) electric dipole moments as functions of temperature. 



coefficient is suppressed by the polaron effect in compar- 
ison with the conventional magnetic Kondo effect. Then, 
Txp decreases rapidly around Tk in sharp contrast to 
Txm around Tk. The shape of Txp at (^i, C/)=(0.2, 1.0) 
around Tk is quite similar to that for (0.2,0.05) [iTj . 
when T is rescaled by Tk. For the case of (0.2,1.0), 
corresponding to the competing region of magnetic and 
electric dipolar Kondo effects, we do not find any sig- 
nificant structure in the entropy and Txm, but in the 
decrease of Txp from the shoulder, the enhanced signal 
can be observed due to the polaron effect in the phonon 
matrix element of xp- 

We remark that the two-channel Kondo effect also 
occurs due to the screening of electric dipole moment 
[H. In fact, at (V"i, C/)=(0. 17384, 0.05), we find a 
plateau of entropy 0.5 log 2 and the significant decrease 
of Txp around the corresponding temperature. Note 
that the shape of Txp at (0.17384,0.05) is different 
from those for (0.2,0.05) and (0.2,1.0). We observe 
the smooth change from Txp=constant to xp=constant 
in this case. Since the electric dipole moment is not 
perfectly screened, vibration still remains and thus, we 
intuitively obtain xp=2/'^2 in the two-channel electric 
dipolar Kondo regime, where 0)2 is another renormalized 
phonon energy, which is different from Cji- It is one of 
future tasks to explain the difference between ui and Cj^ 
by overcoming the difficulty to estimate wi with high 
precision in the NRG calculation. 

Here let us go back to Fig. 1(a). We have found that 
the electric dipole moment is screened in the region of 
large Vi and small U - In order to visualize the change 
of the screened moment, we depict q as the color grada- 
tion in Fig. 1(a). For large Vi and small C/, we find 
the red region with large g, in which electric dipolar 
Kondo effect occurs. On the other hand, there occurs 
magnetic Kondo effect in the blue region of small Vi 
and large U - As a guide of the boundary between mag- 
netic and non-magnetic Kondo regions, we plot inflection 
points of (Pq/dU'^=0 (dotted) and d^q/dV^=Q (broken) 
in Fig. 1(a). For large Vi, two curves run in the green 
area between blue and red regions. Another broken curve 






(b) J 




. 




0.8 




V 




P 




1 




SO.6 


- 


/ 


1 / 


|o.4 

>- 


-^ 


^ 


y 


0.2 


- 


\/ 


/i 


0.0. . — ' L. 


.^ 


. I* . 1 



0.14 0.15 0.16 0.17 0.18 0.19 0.2 0.14 0.15 0.16 0.17 0.18 0.19 0.2 
V, V, 



FIG. 4: (Color online) (a) Sommerfeld constant vs. Vi for 
H—Q and (b) the ratio of the Sommerfeld constants vs. Vi . 



appears on the two-channel Kondo hne for [/<0.6, ex- 
cept for U~0, suggesting the electric dipolar two-channel 
Kondo region. Note, however, that the inflection points 
seem to lose the meaning of boundary in the s-channel 
Kondo region, since renormalized Fermi chain is realized 
in the region of small U and small Vi [12|,[l3| . In the small 
area of yellow and orange of the s-channel Kondo region, 
there also occurs non-magnetic Kondo effect, which is 
interpreted as the Yu- Anderson Kondo effect [lOj. 

Now let us discuss the magnetically robust heavy elec- 
tron state by the Sommerfeld constant 7. For the pur- 
pose, we add the Zeeman term Hz=gsfJ'BH{n-^ — ni)/2, 
where H is an applied magnetic field, to the model (1). 
Then, we evaluate 7 at T^A^^"', the lowest tempera- 
ture at which we can arrive in the present NRG calcu- 
lations. In Fig. 4(a), we show the results for 7 in the 
unit of mJ/mol • K^ with D^l cV. We find that 7 shows 
divergent behavior around the two-channel Kondo fixed 
points, shown by the solid curve in Fig. 1(a). On the 
line of [/=0, we do not find the phase conversion, but 
for FiJ^O.17, the two-channel Kondo line exists in the 
extreme vicinity of U=0. Thus, 7 becomes very large 
for U^O and Vi^^Cl?. In the vicinity of the two-channel 
Kondo line, we find the enhancement of 7 due to the non- 
Fermi liquid properties, irrespective of U. For the param- 
eters away from the two-channel Kondo fixed points for 
C/>0.05, 7 in the region of the magnetic Kondo effect 
is relatively large in comparison with that of the non- 
magnetic Kondo effect of electric dipolc moment. 

In Fig. 4(b), we show the ratio of 7's at H^Q T and 
30 T. If this ratio is near the unity, we judge that 7 is 
magnetically robust. As we easily imagine, in the Kondo 
effect concerning the electric dipole moment, 7 does not 
depend sensitively on the magnetic field. As for criteria of 
the magnetically robust heavy electron state, we consider 
the conditions of 7>1000 and 7(30T)/7(0T)>0.9, shown 
by horizontal broken lines in Figs. 4(a) and (b). After the 
NRG calculations, we find that the magnetically robust 
heavy electron state appears in the region of Vi^O.17 and 
t/^0.1, except for the narrow region near the two-channel 
Kondo line. Note that the ratio is strongly suppressed 



due to the divergent behavior of 7(0T). That region is 
included in the non-magnetic Kondo effect region with 
orange and red in Fig. 1(a). In comparison with the 
mechanism of magnetically robust large 7 on the basis 
of the charge Kondo effect [19| , it seems to be easier to 
obtain large 7 in the present scenario based on the non- 
magnetic two-channel Kondo effect. 

In summary, we have analyzed the two-channel con- 
duction electron model with vibrating magnetic ion. We 
have found two types of Kondo effects due to the alterna- 
tive screening of magnetic and electric dipolc moments. 
Near but not exactly on the two-channel Kondo line with 
electric dipolar origin, we have found that 7 is magnet- 
ically robust. Non- magnetic electric dipolar Kondo be- 
havior is expected to be observed in cage-structure mate- 
rials. In particular, magnetically robust non- Fermi liquid 
behavior is an interesting possibility. 

This work was supported by KAKENHI (20102008). 
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